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(U)  This  Semiannual  Progress  Report  by  Phil^o-^ord ,  Space  and  Re-entry 
Systems  Division,  Newport  Beach,  California,  on  the  Advanced  Penetration 
Problems  Task  of  Wake  Structure  Measurements  was  prepared  for  the  Advanced 
Rt  search  Projects  Agency,  Department  of  Defense,  ARPA  Order  No,  888.  The 
research  was  monitored  by  the  Space  and  Missile  Systems  Organization,  USAF 
under  Contract  F04701-68-C-0032 . 

(U)  The  author  is  Dr,  A.  Demetriades  of  Philco-Ford,  Space  and  Re-entry 
Systems  Division.  The  data  herein  covers  the  reporting  period  of  16  Octo¬ 
ber  1967  through  15  April  1968. 

(U)  This  technical  report  has  been  reviewed  and  is  approved. 

W.  D.  McComb 

Advanced  Penetration  Problems 
Project  Officer 


ABSTRACT 


(U)  Progress  on  the  subject  contract  in  the  period  October  15,  1967  to 
April  15,  1968  is  described  in  this  report.  Work  has  been  performed  on  the 
turbulent  structure  of  a  two-dimensional  supersonic  wake  at  angle  of  attack, 
on  a  similar  wake  at  zero  incidence  but  with  heat  transfer,  and  a  turbulent 
plasma  jet.  Data  from  a  zero- incidence ,  adiabatic  wake  are  being  further 
reduced.  Notable  findings  in  these  experiments  include  highly  accurate 
values  of  the  turbulent  Reynolds  and  Prandtl  number,  a  re-appe  irance  ol  the 
so-called  front  wavelength  of  free  turbulence,  and  the  explanation  cf  ve**y 
large  electron  density  fluctuations  as  due  to  transitional  mixing. 
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SECTION  I 
SUMMARY  OF  WORK 


Reduction  of  data  has  continued  on  the  two-dimensional  turbulent  wake  at  zero 
incidence  and  heat-transfer.  Evaluation  of  the  wake  drag  and  energy  integrals 
along  the  stream  direction  showed  that  these  were  constant  with  distance, 
meaning  that  the  experiment  had  been  "tailored"  properly.  The  virtual  diameter 
was  such  that  a  wake  "map"  to  almost  2,000  diameters  became  possible.  Simi¬ 
larity  in  the  radial  direction  was  demonstrated  by  a  very  consistent  collapse 
of  the  data  on  the  same  curve.  The  expected  agreement  with  low-speed  data  in 
the  far  wake  was  shown,  and,  in  fact,  the  turbulent  Reynolds  number  was  found 
to  be  13.0.  The  matter  of  this  number,  important  to  the  hypersonic  wake, 
should  be  considered  fairly  settled  by  this  time.  In  the  turbulence  aspect, 
the  presence  of  tve  "front  wavelength"  was  again  detected,  bearing  apparently 
the  same  geometric  similarity  to  the  wake  diameter  as  was  found  in  the  axi- 
symmetric  wake. 

The  effects  of  angle  of  attack  and  heat-transfer  on  the  wake  structure,  long 
suspected  to  be  decisive,  have  been  for  the  first  time  tackled  in  compre¬ 
hensive  laboratory  studies.  These  are  performed  with  two-dimensional  geometries 
as  a  logical  first  step  to  introducing  gradually  the  complexities  of  the  final 
problem.  The  experiments  are  now  nearly  complete,  although  the  data  have  not 
been  reduced.  Nevertheless  at  least  one  important  prediction  of  the  Lees-Gold 
stability  theory,  the  delay  of  transition  when  the  wake  is  stabilized  by 
heating,  has  been  verified. 

Further  information  has  emerged  from  the  data  reduction  sequence  of  the  turbu¬ 
lent  plasma  jet  experiment.  It  is  now  known  that  the  front  (interfacial) 
structure  of  the  jet  has  the  same  geometric  features  for  an  electron-sensitive 
observer  that  it  does,  say,  for  a  temperature-sensitive  observer.  The  longi¬ 
tudinal  wavelength  by  which  this  structure  is  organized  is,  apparently,  obeying 
the  same  relation  to  the  jet  diameter  that  the  wake-front  scale  bears  to  the 
wake  diameter.  In  the  zone  of  jet  transition-to-turbulence  the  mixing  process 
has  been  found  to  resemble,  to  some  extent,  the  speculative  models  advanced 
early  in  this  decade  by  Lin  and  Proudian.  Once  the  transition  process  is  com¬ 
plete,  and  for  recombination  reaction  rates  which  are  apparently  too  slow  for 
the  turbulence,  the  electron  density  fluctuation  approach  numerically  the 
fluctuations  in  the  gas  density.  The  integral  scales  of  electron  density  fluc¬ 
tuations  do  not  appear  to  scale  directly  with  the  jet  diameter  as  expected  for 
the  scales  of  neutral  gas  variables.  This  significant  point  is  being  pursued 
further . 

During  the  past  six  months  an  effort  has  been  made  to  review  the  present  ex¬ 
periments  critically  in  the  light  of  work  done  in  other  laboratories,  to 
document  the  present  findings  into  the  form  of  technical  reports,  and  to  dis¬ 
seminate  unclassified  portions  to  the  technical  community  at  large  in  the  hope 
of  eliciting  critical  comments.  Technical  exchanges  have  been  held  with  the 
Caltech  group  under  Prof.  L.  Lees  aad  the  group  at  AC  Electronics  working  with 
Dr.  K.  S.  Wen.  In  the  course  of  technical  presentations  by  RCA/Canada  at 
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SAMSO,  Philco-T'ord  personnel  had  the  opportunity  of  evaluating  experiments 
related  to  the  Philco-Ford  work.  In  April  1968  A.  Demetrades  visited  CAR! 
in  Quebec,  Canada,  and  examined  at  close  quarters  the  work  done  in  the  bal¬ 
listic  range  at  speeds  of  15,000  fps.  A  technical  report  was  written  in 
February  1968  on  the  mean  flow  properties  of  the  turbulent  jet;  two  earlier 
paper  ?  dealing  with  basic  aspects  of  wake  structure  have  been  accepted  for 
publication  in  the  Physics  of  Fluids  Journal  and  the  Journal  of  Fluid  Mechanics, 
and  a  third  paper  dealing  with  the  intermittency  properties  of  the  turbulent 
jet  has  been  approved  for  delivery  at  the  Fluids  and  Plasmas  Conference  of  the 
AIAA  in  June  1968  in  Los  Angeles,  Calif.  On  April  25,  1968,  A.  Demetriades 
delivered  a  talk  entitled  "Structure  of  Compressible  Turbulent  Wakes"  at 
the  12th  AMRAC  meeting  in  Washington,  D.  C.  This  paper  will  appear  in  the 
AMRAC  Proceedings  which  will  be  available  shortly. 


SECTION  II 


STRUCTURE  OF  THE  TWO-DIMENSIONAL  WAKE 


2 . 1  INTRODUCTION 

The  two-dimensional  wake  work  (WED  experiment)  completed  during  the 
previous  contract  period  is  yielding  a  wealth  of  important  information 
through  a  continued  program  of  data  reduction  and  interpretation.  In 
Reference  1  detailed  information  on  this  experiment  is  supplied  os 
regards  the  facilities,  model,  instrumentation  and  procedures,  and  the 
preliminary  evaluation  of  the  mean  and  turbulent  flow  row  data.  In  the 
present  Section  the  final  results  as  regards  the  mean  and  intermittent 
flow  will  be  outlined;  the  turbulence  results  will  be  described  at  the 
present  reduction  stage. 

2.2  MEAN  FLOW  IN  THE  TWO-DIMENSIONAL  WAKE 


2.2.1  ANALYSIS 


A  similarity  analysis  of  the  two-dimensional  compressible  turbulent  wake 
has  been  carried  out  along  lines  similar  to  those  for  the  ax i- symmetric 
wake  (Reference  2).  As  in  all  similarity  analysis  of  comparable  flows, 
the  analysis  utilizes  at  least  one  empirical  fact,  namely  that  the  axis 
velocity  decays  as  the  inverse-square  of  axial  distance.  With  this  input, 
the  results  are  that  the  axis  temperature  defect  also  decays  as  l/>/x, 
and  that  the  density  defect  decay  approaches  the  l/>/x  behavior  much  later 
than  either  the  velocity  or  the  temperature.  The  transverse  scale  I. 
increases  as  Vx,  and  the  radial  profiles  o£  velocity  and  temperature  are 
similar  in  the  radial  parameter  Y/L  where  Y  is  the  Howarth-Dorodnitzyn 
radius . 


Quite  recently  Laufer  (Reference  3)  has  given  a  formal  procedure  ::or 
deriving  the  characteristics  of  mean  compressible  turbulent  (free)  flows 
once  the  corresponding  incompressible  flow  is  known.*  His  results  are 
nearly  identical  to  those  outlined  in  the  previous  paragraph  and  give  the 
following  relation  for  the  axis  velocity  defect  far  in  the  waker, 
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where  u  and  u(o)  are  (in  the  wind-tunnel  coordinate)  the  velocity  just 
external  on  the  wake  and  on  the  axis  respectively,  Rx  is  the  Turbulent 
Reynolds  Number  and  where  x  accounts  for  the  virtual  origin  xc  and  the 
two-dimensional  drag  Cj)h: 
x  -x 

_  _  o 


*D.  J.  Laufer,  head  of  the  Department  of  Aerospace  Engineering  at  the 
University  of  Southern  California,  is  a  Philco-Ford  Consultant. 
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It  is  furthermore  predicted  that  the  factor  in  the  Gaussian  exponent  of 
the  velocity  distribution  is  0.5,  that  is 
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where  r)  =  Y/L. 
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It  should  be  noted  that  the  similarity  law  w  /v  x  differs  in  the^^ 

exponent  in  the  corresponding  law  for  the  axi-symmetric  wake,  w  ~  x 
The  target  of  this  subtask  is  therefore  primarily  the  existence  of 
similarity  and  especially  the  value  of  Rx  which  is  equally  important  to  the 
axi-symmetric  wake  as  well  and  which  is  expected  to  be  insensitive  to 
geometry. 

2.2.2  RESULTS 

Reference  1  presented  experimental  results  from  WED  as  already  reduced 
thru  the  WEB-I  (DR  25F)  program.  In  that  program  quantities  at  various 
radial  positions  were  outputed  as  a  function  of  radius  divided  by  the 
so-called  half-radius  of  the  wake.  This  type  of  presentation  is  a  rather 
insufficient  test  of  similarity,  s*nce  it  forces  all  points  to  pass 
through  too  many  points  -  it  forces,  that  is,  d^ta  at  various  stations  to 
collapse  on  a  single  curve  regardless  of  similarity.  A  new  program, 
WEB-VIII  (SRS  No.  00 6 A)  ,  was  therefore  composed  wh<-se  main  function  was  to 
re-normaiize  the  radius  so  that  quantities  could  be  plotted  versus 
7)  5  Y/L.  The  new  program  also  outputs  the  normalized  physical  radius 
Y/L,  the  integrals  of  momentum  (i.e.  Coh)  and  energy,  the  average  value 
of  Cj)h  and  the  quantity  x,  the  wake  Reynolds  Number  Rew  and  so  on. 

The  energy  integral  and  the  momentum  integral  should  be  constant  along  a 
wake  which  is  unaffected 'by  drag  "swallowing",  that  is  for  a  wake  which, 
as  the  present  experiment  intended,  possesses  only  an  insignificant  shock 
signature.  Figure  1  shows  that  indeed  this  is  the  cose.  The  average 
Cj)h,  computed  by  excluding  the  first  few  points  which  are  obviously 
affected  pc  tly  by  instrument  resolution  is  0.00909  cm.  This  length  can 
be  thought  of  as  the  effective  frontal  height  of  a  two-dimensional  body  of 
drag  coefficient  of  unity.  It  is  not  coincidental  perhaps  that  the 
actual  frontal  height  of  the  ribbon  is  0.01  cm,  within  10  percent  of 
chat  figure. 

The  radial  distributions  of  velocity  and  temperature  are  shown  on 
Figure  2  as  plotted  versus  Tj  .  The  correlation  is  quite  satisfactory. 

It  is  especially  significant  that  the  velocity  data  u  pass  through  the 
point  (u  =  0.5,  T)  =  1)  quite  unaided,  signifying  a  near-equality  of  the 
scale  L  and  the  half-radius.  The  simi1arity  distribution  for  u  and  T: 
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FIGURE  I  VARIATION  OF  TWO-DIMENSIONAL  WAKE  MOMENTUM  AND  ENERGY  THICKNESS 
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FIGURE  2  RADIAL  VARIATION’S  OF  NON-DIMENSIONAL  WAKE  VELOCITY  (TOP) 

AND  AMBIENT  TEMPERATURE  (BELOW) 
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T  =  T(o)  -  T„ 
is  Gaussian  in  character. 


(4) 


The  axial  variations  of  the  scale  L  and  of  the  defects  have  already  been 
presented  in  Figures  34,  35  and  36  of  Reference  1.  Since  the  abscissa 
XSTATION  of  these  figures  ij  simply  proportional  to  the  correct  axial 
parameter  x,  it  is  indeed  see”  that  L  'vy/lf  and  w  ,0  ~  x  as  expected. 

We  can  now  present, in  addition,  the  value  of  the  Reynolds  Number  R?  as 
obtained  from  the  present  experiment.  Using  equation  (1): 


■*T 
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we  obtain  the  results  show  on  Figure  3.  The  average  value  (again 
excluding  a  short  distance  near  the  body)  is  13.0,  remarkably  close  to  the 
12.7  found  by  Tounsend  at  low  speeds  (Reference  4)  and  to  the  12.8  found 
for  the  supersonic  aci- symmetric  wake  (Reference  2). 

The  identity  among  these  numbers  is  of  course  hardly  surprising  in  view  of 
the  fact  that  all  wake  flows  eventually  decelerate  to  the  speed  of  the 
ambient  environment.  Intense  compressibility  effects  cannot  be  extracted 
from  the  present  experiment.  In  fact  spot-checks  of  the  effect  of  plotting 
the  radial  variations  versus  the  physical  radius  Y/L  (rather  than  the 
Y/L)  were  made.  Typically  the  scatter  in  Y/L  at  u  =  0.5  was  about 
15  percent  of  the  scatter  in  ^/L  -  too  small  an  effect  from  which  to 
draw  conclusions. 

2.3  INTERMITTENCY 

The  intermittent  properties  of  the  two-dimensional  wake  have  been  measured 
exactly  by  the  same  procedure  as  used  for  the  WEB  (see  Reference  5) .  The 
intermittency  factory  has  been  found  to  be  distributed  normally  along  the 
wake  radius,  i.e.,  y  is  a  Gaussian  function  of  radial  distance  about  the 
position"Y  of  the  turbulent  front.  Figure  4  shows  the  growth  of”  and  of 
the  standard  deviation  a  with  distance;  as  expected,  Y  and  a  are 
approximately  growing  as  the  square  root  of  distance.  This  growth  is 
shown  better  in  Figure  5,  where  these  two  lengths  have  been  normalized  by 
the  ‘ransverse  scale  L.  The  average  value  of  Y/L  is  about  1.96,  which  is 
only  7  percent  than  the  value  obtained  for  the  ax i- symmetric  wake. 

However  the  standard  deviation  has  an  average  value  of  about  0.38L,  which 
is  considerably  lower  (by  a  factor  of  2)  than  the  standard  deviation  of 
the  axi- symmetric  wake. 
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FIGURE  5  NON-DIMENS 


2.4  THE  TURBULENT  FLOW 


In  Reference  1  the  description  of  the  turbulence  flow  results  obtained 
with  the  WEB  experiment  was  limited  to  the  raw  data.  At  this  time  the 
reduction  of  the  latter  data  has  progressed  through  the  WEB-II  and  parcly 
through  the  WEB-III  programs. 

The  objective  of  the  WEB-II  program  (Reference  6)  is  to  calculate  the 
response  of  the  overall  electromechanical  system  (including  the  hot  wire) 
at  each  frequency  passband  and  to  thus  produce  a  numerical  factor  by  which 
the  raw  signal  in  that  band  is  to  be  corrected.  This  program  then  outputs 
the  correct  fluctuation  spectrum  of  the  wire  output,  although  further 
resolution  of  the  spectrum  into  velocity  and  density  contributions  (the 
so-called  Modal  resolution)  occurs  later.  On  the  basis  of  this  spectrum, 
three  relevant  items  will  be  discussed  in  this  section:  the  extent  of 
turbulence  energy  along  the  spectrum,  the  differences  between  uncorrected 
and  corrected  spectra  and  the  appearance  of  the  spectral  prominence 
thought  to  represent  the  front  wavelength. 


The  Reynolds  number  of  the  wake  is  defined  by 
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where  V#  is  the  kinematic  viscosity  of  the  receiving  medium.  In  the 
two-dimensional  wake  this  number  is  constant  along  the  wake  because  the 
velocity  difference  Uqo  -  u(o)  decreases  as  while  L  increases  as  /*. 
In  the  present  experiment  Rew  =  158.  Now  there  is  another  relevant 
Reynolds  Number,  that  of  turbulence 


Re 


t 


(7) 


where  Au  is  rms  local  velocity  fluctuation,  A  the  integral  fluctuation 
scale  and  V  the  mean  local  kinematic  viscosity.  From  the  WEA  -  WEB 
experiments  (Reference  2)  it  is  known  that  A  a  L  and  that  Au  a  0.4 
(Uqj,  -  u(o)j  on  the  wake  axis;  far  in  the  wake  va v^.  As  a  result  the 
following  estimate  of  Refc  can  be  made: 


Re 


t 


Au.  A 
v 


=  0.4  Re  a  63 
w 


(8) 


Furthermore,  the  stream  velocity  in  the  far  wake  is  of  the  order 

60,000  cm/sec  and,  from  the  WEB-VIII  program,  that  L  is  of  the  order 

0.08  cm  in  the  far  wake.  Thus  the  frequencies  registered  by  the  hot-wire, 
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as  eddies  of  integral  scale  size  pass  it,  are  of  the  order  60,000/0.08, 
that  is  750  kilohertz.  In  turbulence  there  exist  eddies  considerably 
smaller  than  the  integral  scale,  so  that  the  spectral  content  being 
measured  should  extend  into  several  megahertz. 

This  highly  extended  frequency  range  is  too  wide  to  be  handled  in  its 
entirety  by  the  present  hot-wire  system.  The  system,  it  should  be  recalled, 
consists  of  state-of- the-art  componentry,  is  aided  by  the  forementioned 
computer  program  (WEB-II)  and  it  can  reliably  return  results  in  the  range, 
say,  to  about  1  MHZ.  The  net  result  is  that  in  the  present  work  no  attempt 
was  made  to  go  beyond  about  1.2  MHZ  in  the  spectra,  and  the  turbulence 
content  for  the  small-eddy  sizes  was  not  measured.  However,  important 
features  were  seen  in  the  range  studied,  and  the  frequency- integrates  but 
modally-resolved  fluctuations  should  be  correctly  deduced  since  not  much 
energy  is  anyway  contained  in  the  small  scales. 

With  these  notations,  consider  the  spectra  shown,  for  example,  on  Figure  6. 
The  objective  here  is  to  consider  the  effect  which  the  restoration  of 
frequency  response  effected  by  the  WEB-II  program  lias  on  the  spectral  shape. 
Observe  that  the  response  inadequacies  of  the  system  without  response 
restoration  decreases  tremendously  the  apparent  energy  content  at  the  higher 
frequencies.  For  example  in  the  vicinity  of  1  MHZ,  the  uncorrected  spectrum 
shown  is  down  100,000  times  the  zero- frequency  value  (in  the  square  of 
spectral  density)  whereas  with  proper  restoration  this  decrease  is  only  by  a 
factor  of  2.  The  response-restoration  technique  thus  is  crucial  to  the 
present  experiment. 

A  final  point  concerns  the  peak  shown  on  Figure  6  and  more  adequately  on 
Figure  7.  This  sort  of  energy  concentration  is,  formally  speaking,  quite 
unexpected  in  the  one-dimensional  spectra  of  turbulence.  Similar  peaks  were 
observed,  however,  in  the  one-dimensional  spectra  of  the  axi-symmetric  wake 
and  have  already  been  described  in  Reference  6.  In  that  work  the  peak  was 
explained  as  arising  from  the  periodic  structure  of  the  turbulent  front  which 
is  superposed  on  the  bona  fide  turbulence  spectrum  in  the  region  of  intense 
intermittent  activity,  i.e.,  away  from  the  axis.  If,  as  shown  in 
Reference  5,  this  signal  has  a  wavelength  related  to  the  wake  diameter  (or 
the  scale  L;  one  is  led  to  believe  that  the  signal  would  appear  in  WED  at  a 
much  higher  frequency,  because  the  WED  diameter  is  much  smaller  than  the  WEA 
diameter.  Indeed,  the  peak  frequencies  shown  on  Figure  6  and  7  occur  at 
about  400  KHZ,  whereas  the  peaks  in  WEA  appeared  at  about  50  KHZ.  This 
ratio  of  8  in  the  wavelength  is  about  equal  to  the  ratio  of  the  axi- 
symmetric  wake  diameter  to  the  two-dimensional  wake  diameter. 


-12- 


B0865  U 


FIGURE  6 


UNRESOLVED  SPECTRA  OF  TURBULENCE  AT  A  FIXED  POINT  IN  THE  WAKE, 
ILLUSTRATING  THE  EFFECT  OF  RESPONSE  RESTORATION (PROPER 
FREQUENCY  RESPONSE) 


SPECTRAL  INTENSITY 
(ARB  UNITS) 


FREQUENCY  (KHz) 

FIGURE  7  RESTORED,  ONE -DIMENSIONAL.  UNRESOLVED  SPECTRA  AT  FIXED 

radial  but  varied  axial  positions 
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SECTION  III 


THE  TURBULENT  PLASMA  JET 


During  the  present  contract  period  the  work  on  the  mean  and  intermittent 
properties  of  the  turbulent  argon  plasma  jet,  including  the  reduction  of 
data,  has  been  completed.  Data  reduction  is  proceeding.  This  section 
describes  in  final  form  the  results  of  the  intermittency  measurement  and 
the  continuing  work  on  the  turbulence  measurements.  The  render  is  referred 
to  reports  from  the  previous  contract  period  for  details  on  the  flow 
facilities,  instrumentation  and  procedures. 

3.1  THE  INTERMITTENT  FLOW 

3.1.1  DEFINITIONS  AND  PROCEDURES 

The  measurement  of  the  intermittent  flow  in  the  jet  is  a  natural  and  yet 
important  step  in  understanding,  eventually,  the  interplay  between  gas- 
dynamics  and  chemistry  in  the  turbulent  plasma.  As  can  well  be  imagined, 
two  important  questions  arise:  (1)  how  does  the  high  temperatures  in  the 
jet  affect  the  statistically  averaged  geometric  features  of  the  front? 

(2)  how  does  the  front  geometry  and  behavior  appear  to  an  observer  sensitive 
to  the  electrons  vis-a-vis  an  observer  sensitive,  say,  to  the  temperature 
of  the  flow? 

In  the  present  experiment  the  front  properties  were  measured  with  the  hot¬ 
wire  anemometer  and  the  Langmuir  probe.  The  signal  from  the  hot-wire 
amplifier  (or  the  ac  amplifier  connected  across  the  load  of  the  Langmuir 
probe)  was  fed  to  the  intermittency  circuit. 

The  intermittency  factor  y  was  recorded  at  each  of  about  20  radial  positions 

Y  at  each  x.  The  rms  front  position  Y  and  its  standard  deviation  o  were 
then  obtained  from  the  usual  formulas: 

V  =  |j  (Y  -  Yr  dCY  *  Y)  1/2 

In  these  formulas  computational  convenience  requires  that  y  is  expressed  as 
a  function  ct  Y.  This  was  done  by  polynominal-fitting  the  experimental 
points  with  .'he  computar. 

In  addition  to  Y  and  O  the  autocorrelation  scale  of  the  front  is  of  some 
interest.  If  the  front  behaves  like  a  stationary  random  variable  it  is 
possible  to  find  this  scale  from  a  measurement  of  the  so-called 


(8) 

(9) 


-15- 


I 

I 


null-crossings  (or  null  frequencies)  of  the  front,  that  is  by  the 
frequency  N0  =  N(Y)  at  which  a  sensor,  located  at  Y  =  Y,  enters  the  front. 

To  this  end  the  distribution  of  null  frequencies  N(Y)  with  radius  were 
recorded  by  electronically  counting  the  rectangular  pulses  at  the 
Schmidt-trigger  output  of  the  intermittency  circuit.  These  distributions, 
which  of  course  peak  at  or  near  Y,  were  also  computer-fitted  by  polynominals. 
The  microscale  of  the  front  is  given  by 


(10) 


where  u(Y)  is  the  jet  velocity  at  Y.  An  alternate,  and  possibly  more  mean¬ 
ingful  front  "wavelength"  is  simply 


(ID 


3.1.2  RESULTS  AND  DISCUSSION 


Figure  8  shows  a  typical  comparison  of  the  distributions,  with  radius,  of 
the  intermittency  factor,  the  temperature,  and  the  electron  density.  As 
expected,  the  fluid  is  strongly  intermittent  at  distances  from  the  axis  at 
which  the  mean  temperature  is  very  close  to  the  ambient.  Also  from 
Figure  8  it  can  be  seen  that  the  corresponding  relation  between  electron 
density  and  the  intermittency  factor  is  analogous  but  highly  inconsistent 
from  one  x  to  another:  at  a  fixed  value  of  7  =  0.5,  for  example,  n  varies 
considerably  in  the  x-range  investigated. 

Important  questions  to  be  settled  in  analyzing  the  intermittency  data 
include  the  randomness  of  the  interface,  the  coincidence  between  electron- 
and  temperature-interface  and  the  comparison  of  the  front  geometry  and 
behavior  to  that  of  low-speed  jets.  The  growth  of  the  front  "radius" 

Y  is  shown  dimensionally  in  Figure  9  and  non-dimensionally  as  Y/L, 
together  with  the  standard  deviation  o/L,  in  Figure  10.  There  is  little 
dout  that  these  two  quantities  grow  lirearly  and  scale  with  L,  and  that 
they  certainly  are  independent  of  the  method  of  measurement.  With  Y 
and  a  so  computed  it  is  possible  to  make  a  test  of  the  randomness  of  the 
front  such  as  is  shown  on  Figure  11;  it  is  seen  that  the  Gaussian 
distribution  of  the  intermittency  factor  about  the  front  position,  that  is 
with  (Y-Y) /a,  again  obtains  in  the  fashion  characteristic  of  stationary 
random  variables.  In  the  same  figure  further  evidence  of  randomness  is 
supplied  by  the  distribution  of  the  normalized  crossing  frequencies 
N/k’max  about  the  radial  location  Ymax  =  Y(Nmax)  .  The  scatter  is 
appreciable,  but  enough  data  points  are  available  to  illustrate  the 
Gaussian  behavior  of  N/Nmax  characteristic  of  randomness.  As  before, 
there  is  no  systematic  difference  between  data  taken  with  the  hot-wire 
and  data  taken  with  the  Langmuir  probe. 
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ELECTRON  DENSITY  DISTRIBUTIONS 
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FIGl’RE  10  THE  TURBULENT  JET  FRONT  POSITION  AND  ITS  STANDARD  DEVIATION 
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DISTRIBUTION  OF  NULL  FREQUENCIES 
ABOUT  THE  FRONT  POSITION 
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FIGURE  11  DISTRIBUTION  OF  THE  INTERM ITTENCY  FACTOR  AND  NULL-CROSSING  FREQUENCIES  ABOUT  THE 
FRONT  POSITION 


-20- 


The  null  frequencies  N0,  defined  as  the  crossing  frequencies  measured  at 
the  position  of  the  interface,  were  found  to  be  very  close  to  the  maximum 
frequency  Nmax  measured  (at  each  x)  at  Ymax  thus,  N0 =  N(Ymax)  =  N(?)  = 

Nmax  as  expected.  Figure  12  shows  the  variation  of  N0  with  x;  again,  the 
frequencies  measured  with  the  hot-wire  and  the  Langmuir  probe  are  about  the 
same  at  each  x.  The  decrease  of  N0  shown  in  this  figure  immediately  suggests 
the  lengthening  of  the  front  scale  \p  (or  Ap)  expected  as  the  turbulent  flow 
widens  with  distance.  However,  the  computation  of  these  scales  from 
Equations  (10)  and  (11)  cannot  be  performed  with  the  velocity  data  shown  so 
far;  these  latter  are  mean  velocities,  whereas  what  is  of  course  needed  is 
the  velocity  of  the  front  itself  or  at  least  the  mean  velocity  within  the 
front  -  the  so-called  eddy  velocity. 

Measurement  of  eddy  velocities  within  "corrugations"  in  the  jet  boundary 
were  made  with  the  double  Langmuir  probe  (two-wire  correlation  probe). 

Eddy  velocities  with  such  a  probe  can  best  be  measured  from  the  cross¬ 
correlation  function  of  the  probe  signals;  in  the  present  work  the  more 
primitive  method  was  adopLed  of  obtaining  a  finite  number  of  dual-trace 
oscillograms  of  the  turbulence  pattern  as  it  was  convected  by  the  wires. 

From  these  samples  a  rudimentary  statistical  distribution  of  the  local  eddy 
velocity  was  constructed  at  each  point  from  which  the  mean  velocity  (and 
its  fluctuation)  could  be  extracted.  The  (mean)  eddy  velocity  at  Y  was 
then  used  in  place  of  u  in  Equations  (10)  and  (11). 

Reference  5  describes  experiments  done  on  the  intermittency  features  of 

compressible  wakes  and  advises  the  use  of  the  "front  wave  length"  Ap  (rather 

than  \f)  as  a  more  significant  measure  of  the  front  geometry.  In  doing  so, 

evidence  is  drawn  from  observations  of  the  turbulence  spectra  (near  the 
front  posicion  Y)  that  the  front  itself  may  not  be  truly  random,  but  rather 
weakly  periodic  at  a  wavelength  A p.  The  asymptotic  fromAp/L  =  u/N  L  of 
this  wavelength  for  the  axi-symmetric  wake  is  dravm  in  Figure  13  in  which 
the  jet  front  wavelength,  obtained  per  the  previous  paragraph,  are  also 
presented.  The  agreement  between  these  two  independent  experiments  is  quite 
revealing.  The  wavelength  of  the  jet  front  appears  to  bear  the  same 
relation  to  L  (about  p/L  s  9)  that  it  does  for  the  wake. 

It  should  be  noted  that  the  physical  dimensions  of  the  front  are  given  here 
in  physical  coordinates,  and  this  should  be  accounted  for  comparing  the 
present  results  with  other  experiments.  In  the  latter  the  radius  is  usua  Hy 
normalized  with  the  so-called  half-radius  *\/2-  Thus  Becker  et.al 
(Reference  7)  gives  Y/r^/2  =  1-78,  Bradbury  (Reference  8)  1.73  and  Corrsin 
and  Kistler  (Reference  9)  values  ranging  from  1.65  to  2.00.  From  the  present 
experiment  (see  Figure  8)  we  obtain  Y/L  =  3.35  on  the  average,  and  in  terms 
of  the  half-radius  we  get  Y/r^/2  =2.0  It  therefore  appears  that  quite  a 
few  jet  experiments  generally  agree  as  to  Y/r^/2-  It  should  be  remembered, 
however,  that  heating  the  fluid  generally  increases  the  jet  radius;  the  end 
result,  shown  on  Figure  14,_is  that  our T  extends  farther  out  by  about 
30  percent,  than  Corrsin' s  Y  for  example. 
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WAVELENGTH  OF  TURBULENT  FRONT 


FIGURE  13  VARIATION  OF  FRONT  WAVELENGTH  OF  TYPICAL  FREE  TURBULENT  FLOWS 
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FIGURE  14  COMPARISON  OF  TURBULENT  FRONT  GROWTH  WITH  TYPICAL  LOW-SPEED  DATA 


The  vfl hie  of  a/Y  obtained  here  is  0.23,  not  far  from  Bradbury's  value  of 
/v  °n  Corrsin  s  ranBC  of  0.23  to  0.29.  However  Becker  et.al.  obtained 
Q/Y  »  0.16.  This  might  be  interpreted  as  a  difference  in  intermittent 
properties  between  an  inert  contaminant  and  the  temperature  field,  but  it  is 
more  likely  that  Becker's  method  of  measurement  and  perhaps  differences  in 
the  turbulent  Reynolds  number  are  responsible.  Becker  and  his  co-workers 
have  at  any  rate  obtained  a  unique  relation  between  the  intermittcncy  factor 
and  the  contaminant  density  distribution  with  radius.  In  the  reacting  jet 
employed  here  the  electrons  are  distributed  radially  in  an  irregular  manner 
because  of  reactions  (see  Figure  8)  and  this  relation  does  not  materialize. 

3.2  THE  LARGE  ELECTRON  DENSITY  FLUCTUATIONS 

Perhaps  the  single  most  interesting  phenomenon  observed  during  the  turbulence 
measurements  was  the  very  large  magnitudes  of  the  electron  density  fluctua¬ 
tion.  It  is  statistically  possible  for  the  standard  deviation  of  a  stationary 
variable  to  exceed  the  mean  value  of  the  variable  when  the  amplitude  distri¬ 
bution  of  the  latter  is  extremely  skew.  Such  situations  exist  in  the 
intermittent  zones  of  ordinary  turbulent  jets  for  any  flow  variable  Q  which 
fluctuates  between  nil  in  the  ambient  ("irrotational")  environment  and  some 
mean  level  within  the  turbulent  front.  If  the  latter  does  not  decrease  too 
fast  with  increasing  radius  then  i.:  is  well  known  thatAQ/Q  can  diverge 
with  radial  distance.  In  fact  Becker,  Hottel  and  Williams  (Reference  7) 
have  observed  this  divergence  by  utilizing  an  inert  contaminant  concentration 
for  an  observable. 

The  possibility  that  the  large  fluctuations  on  the  axis  pictured  on 
Figure  15  are  caused  by  the  intermittcncy  is  excluded-by  the  simple  fact 
that  the  flow  on  the  axis  is  fully  turbulent  as  already  demonstrated. 

From  the  results  thus  far  shown  in  fact,  it  appears  that  we  have  here  a  fully 
turbulent  fluid  with  normal  temperature  fluctuation  but  with  highly-skewed 
electron  fluctuations.  An  osci 1 loscopic  study  of  the  turbulent  jet 
immediately  showed  large  "spikes"  in  the  electron  density  time  history  at  a 
point;  these  spikes  are  indeed  what  statistics  would  require  to  produce  the 
large  An/n  of  Figure  15.  Figure  16  shows  a  comparative  oscilloscope  study 
of  the  Langmuir  probe  and  hot-wire  output  at  various  radial  positions  but 
at  the  same  oistance  from  the  jet  nozzle.  By  properly  adjusting  the  scope 
gain  it  is  possible  to  show  conclusively,  first,  the  differences  between 
intermittent  and  fully  turbulent  flow  a^d,  secondly,  the  differences  between 
Langmuir-probe  ai.d  hot-wire  output  in  the  fully  turbulent  region.  Attention 
is  drawn  to  the  large  spikes  in  the  Langmuir  probe  traces  which  are  clearly 
absent  from  the  hot-wire  traces. 


The  appearance  of  the  large  spikes  persisted  for  several  inches  in  the 
upstream  end  of  the  turbulent  zone  studied,  but  decayed  and  eventually 
vanished  as  the  probe  was  moved  downstream;  the  fluctuation  level An/n  also 
decayed  below  unity  as  Figure  15  shows.  To  complete  the  picture,  an 
investigation  of  the  laminar  and  transitional  portions  of  the  jet  was  made 
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FIGURE  15  VARIATION  OF  NORMALIZED  ELECTRON  DENSITY  FLUCTUATIONS  ALONG  THE  JET  AXIS 
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FIGURE  H  OSCILLOSCOPIC  STUDY  OF  TEMPERATURE  AND  ELECTRON  DENSITY  FLUCTUATIONS  (1  OF  2) 
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FIGURE  16  OSCILLOSCOPE  STUDY  OF  TEMPERATURE  AND  ELECTRON  DENSITY  FLUCTUATIONS  (2  JF  2) 


using  a  Langmuir  probe  thicker  than  described  previously  in  this  paper,  so 
that  it  could  withstand  the  much  higher  heat  fluxes  in  the  laminar  jet.  For 
obvious  reasons  the  output  of  this  probe  was  not  quantitatively  interpretable, 
but  it  made  possible  to  piece  together  a  reasonable  description  of  the  events 
leading  to  the  generation  of  the  large  fluctuations. 

Figure  17  shows  this  picture.  The  Langmuir  probe  observation  made  is 
typified  by  an  instantaneous  "snapshot"  of  the  electron  density  on  the  jet 
axis  from  the  nozzle  exit  well  into  the  turbulent  zone.  In  the  laminar  jet 
the  mean  electron  density  is  high  with  an  unsteady  component  consisting  of 
the  rectifier-driven  360  HZ  modulation  whose  rms  magnitude  is  about  10 
percent  of  the  mean.  At  the  flow  velocities  prevalent  here  the  wavelength 
of  this  sinusoidal  fluctuation  is  actually  very  much  larger  than  the  length 
of  the  laminar  jet  (about  10  feet  compared  to  about  1  foot).  In  the  transi¬ 
tion  region  the  breakdown  process  mixes  ambient  fluid  with  the  ionized  fluid; 
the  former  is  electron-poor  so  that  large  "negative"  pulses  (if  one  thinks  of 
the  electron  current  as  producing  a  positive  displacement  of  the  scope 
beam)  appear  as  shown  on  Figure  17:  the  signal  is  thus  "negatively 
cne-sided".  The  mean  electron  density  thus  decreases  below  the  maximum 
level  -  which  also  decreases  because  of  recombinations.  Farther  downstream 
the  spikes  widen  (but  do  not  necessarily  deepen),  and  still  farther  the 
signal  begins  reversing  itself,  i.c.,  becoming  "positively  one-sided". 

In  this  latter  case  the  mean  electron  density  drops  considerably  below  the 
tops  of  the  spikes,  the  fluctuation  distribution  is  very  skew  and  the  rms 
fluctuation  level  is  considerably  higher  than  the  mean  level.  It  is  at  this 
stage  that  the  oscillograms  showi  on  Figure  16  were  taken. 

It  should  be  again  emphasized  that  at  this  stage  the  fluid  is  fully 
turbulent  from  the  standpoint  of  intermi t tency ,  and  that  turbulent  fluctua¬ 
tions  in  the  electron  density,  of  orde^  2o  to  40  percent  of  the  mean,  exist 
in  addition  to  the  large  electron  densi  pikes.  The  latter  are  obviously 
the  remnants  of  the  laminar  flow  upstrear,.  of  the  transition  zone,  in  which 
the  electron  density  level  is  very  high  when  compared  to  the  level  in  the 
turbulent  fluid  itself.  Eventually  these  spikes  also  disappear  and,  as 
Figure  15  shows,  the  ratio  An/n  settles  to  values  of  order  20  percent  in  the 
far  jet. 

In  Figure  17  it  is  schematically  shown  that  the  maximum  of  the  normalized 
fluctuations  is  attained  where  the  mean  electron  level  is  much  below  the 
maximum  level  n^,  that  is  the  locus  of  the  spike  maxima.  If  the  spikes 
are  truly  the  results  of  the  procr ss  of  engulfment  of  the  ambient, 
electron-poor  fluid  by  the  electron-rich  fluid  of  the  jet,  and  if  the 
electron  density  of  the  former  is  nc ,  then 
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FIGURE  17  RECONSTRUCTED  "SNAPSHOT"  OF  THE  DISTRIBUTION  OF 
ELECTRON  DENSITY  ALONG  THE  JET 
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where  P/m  is  the  ratio  of  the  total  mass  to  the  jet  mass.  Here  it  has 
been  assumed  that  both  fluid  components  are  laminar  so  that  the  signal 
can  locally  oscillate  between  nm  and  nc  only,  in  the  manner  suggested 
some  time  ago  by  Feldman  and  Proudian  (Reference  10)  and  Lin  (Reference  11). 
The  type  of  fluctuations  of  Equation  (12)  are  nil  before  entrainment  begins 
(P/m  =  1)  and  also  when  P/m  tends  to  infinity  in  the  far  jet,  and  they 
peak  in  between  with  a  magnitude  depending  on  nm/nc  For  equilibrium 
chemistry,  where  and  nc  are  directly  obtainable  from  the  corresponding 
temperature  Tm  and  Tc ,  it  is  shown  by  Reference  12  that  An/n  of  order  10° 
or  higher  can  obtain  for  practical  values  of  Tm/Tc.  For  finite  reaction 
rates  the  ratio  An/n  becomes  extremely  sensitive  to  the  chemistry. 

Nevertheless  the  simplicity  of  the  jet  flow  brought  about  by  gas-dynamical 
self-preservation  allows  certain  broad  tests  of  Equation  (12)  to  be  made 
easily.  For  example  the  entrainment  rate  is  given  theoretically  and  has  also 
been  measured  directly.  For  recombination  reactions  the  nm  might  also  be 
obtained  as  a  function  of  x.  The  normalized  electron  density  fluctuations 
of  Equation  (12)  thus  become  dependent  only  on  the  ratio  nm0  (of  the 
maximum  electron  density  in  the  laminar  flow)  to  the "ambient"  value  nc. 

The  former  is  actually  obtainable  from  the  experiment  but  the  latter  is 
not,  since  the  electron  density  at  the  base  of  the  spikes  was  detected 
to  be  turbulent  and  of  a  mean  value  much  higner  than  that  outside  the  jet. 
Computations  made  with  Equation  (12)  using  estimated  values  of  nc  gave 
reasonable  numerical  agreement  with  the  measured  axial  variation  of 
An/n  shown  on  Figure  15. 

Equation  (12)  thus  appears  to  be  a  reasonable  descriptor  of  the  physical 
mechanism  giving  rise  to  the  large  electron  density  fluctuations  -  on  the 
strength  of  the  mixing  model  it  describes.  It  does  not  necessarily  endorse 
the  "laminar  mixing"  ideas  of  References  10  and  11  since  the  participating 
lumps  of  fluid  (whether  electron-rich  or  electron-poor)  are  turbulent 
in  character.  At  this  time  it  also  appears  to  affect  the  transitional 
zone  more  so  than  the  turbulent  activity  in  the  far  jet;  in  the  present 
experiment  thus  zone  is  only  a  few  jet  thicknesses  long. 


» 
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SECTION  IV 


SUPERSONIC  WAKES  AT  ANGLE  OF  ATTACK 


4.1  INTRODUCTION 

Drastic  changes  are  introduced  in  the  wake  of  a  body  flying  at  supersonic 
speeJs  when  the  body  is  inclined  at  soir°  incidence  (angle  of  attack)  with 
respect  to  the  flight  vector:  first,  the  increased  total  drag  generates  an 
intense  inviscid  wake,  i.e.,  large  gradients  appear  in  the  flow  surrounding 
the  inner,  or  turbulent  wake.  Secondly,  these  gradients  are  generally  asym¬ 
metric  about  the  wake  axis.  Thirdly,  changes  in  the  viscous  drag  change  the 
characteristics  of  the  viscous  wake.  Finally,  with  the  generation  of  lifting 
forces  vortical  flows  migh  also  appear,  causing  added  complexities  to  the 
flow  field. 

Although  the  new  character  of  the  turbulent  wake  is  hard  to  assess  in  these 
circumstances,  its  role  in  the  mixing  process  can  be  predicted  along  general 
lines.  Its  main  complexity  should  be  its  growth  in  and  interaction  with  the 
asymmetric  inviscid  wake.  The  flow  thus  consists  of  a  "shear  layer"  contri¬ 
buted  to  by  the  inviscid  wake,  with  a  superposed  wake  flow.  It  is  quite 
clear  that  the  behavior  of  this  composite  flow  cannot  always  follow  simple 
wake  rules  of  lateral  growth,  for  example,  or  of  velocity  decay,  since  shear 
layers  and  wakes  are  generally  governed  by  different  rules.  Thus,  a  pure 
shear  layer  supplies  a  constant  reference  velocity  scale  by  definition,  while 
the  corresponding  velocity  scale  of  the  wake  decreases  with  distance.  In 
limiting  cases,  the  competition  between  the  shear-layer-like  and  wake-like 
behavior  can  be  foreseen:  for  small  angles  of  attack  and  large  shear  drag, 
the  wake-like  behavior  should  dominate  early,  while  for  inclination  angles 
large  enough  to  cause  very  persistent  inviscid  asymmetries  the  flov;  is 
probably  shear-layer-like  for  large  distances. 

The  present  experiment  was  designed  mainly  in  a  way  so  as  to  accentuate  this 
competition  between  the  inviscid  and  the  viscous  (turbulent)  wake.  Calcu¬ 
lations  of  the  supersonic  inviscid  flow  behind  a  two-dimensional  slender 
body  at  incidence,  verified  experimentally,  showed  a  considerable  persistence 
of  asymmetries  in  a  pattern  approaching  an  ideal  shear  layer.  Thus,  the 
problem  is  not  only  relevant,  but  was  done  in  a  way  capable  of  fundamental 
conclusions  on  the  basic  behavior  of  fluids. 

4.2  MODEL  AND  FACILITIES 

This  vjork  is  performed  in  the  Philco-Ford  supersonic  wind-tunnel  at  a  con¬ 
tinuous  air  flow  at  Mach  3.0,  a  stagnation  pressure  of  730  mm  Hg  absolute 
and  a  total  temperature  of  39°C;  the  tunnel  is  described  in  detail  in 
Reference  13.  The  model  consisted  of  a  0.004-inch  thick,  0.116-inch  wide 
stainless-steel  ribbon  stretched  across  the  tunnel  test-section,  identical 
to  that  used  in  the  WED  experiment  (two-dimensional  adiabatic  wake  at  zero 
incidence)  and  described  in  detail  in  Reference  1.  A  major  modification  to 
the  model  support  consisted  of  the  addition  of  a  fixed  steel  protractor  and 
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pointer  by  which  incidence  angle  could  be  set  and  resolved  to  better  than 
0.25  degree  twist-angle  resolution  is  about  0.1  degree.  Other  modifica¬ 
tions  to  the  support,  described  in  Section  V  were  made  to  accomodate  the 
heat-transfer  experiment.  The  protractor  is  shown  on  figure  18.  figure  19 
shows  the  arrangement  of  the  model  and  the  nomenclature  of  the  flow  phenomena 
and  coordinate  axes. 

A. 3  EXPERIMENT  DESIGN 


A. 3.1  INVISCID  WAKE 

Detailed  calculations  of  the  inviscid  wake  of  the  flat  plate  were  mad<*  using 
the  method  of  two-dimensional  characteristics  at  an  angle  of  attack  Oi 
20  degrees;  for  different  angles,  the  flow  field  should  not  alter  greatly 
from  this  calculation,  provided  that  the  incidence  angle  for  shock  detach¬ 
ment  at  the  leading  edge  (about  3A  degrees  at  this  Mach  number)  is  not 
exceeded.  In  this  calculation  the  upstream  boundary  conditions  on  the  wind¬ 
ward  side  were  the  oblique  shock  at  the  leading  edge;  for  the  leeward  side, 
the  trailing  edge  shock  as  modified  by  the  leading— edge  expansion  fan. 

Figure  20  shows  graphically  the  results.  Attention  should  be  drawn  to  tue 
two  nearly  parallel  streamlines  showning  on  either  side  of  the  wake  axis. 

It  turns  out  that  these  demark  abruptly  changing  conditions  in  the  sense, 
of  flow  gradients  (not  in  the  sense  of  flow  discontinuities).  Ihe  net  re¬ 
sult  is  that  between  each  of  these  two  streamlines  and  the  axis  there  appears 
a  uniform  "plateau"  in  flow  conditions;  the  two  adjacent  "plateau"  constitute, 
then,  a  shear-layer  situation  which  apparently  does  not  decay  much,  either  in 
extent  or  intensity,  as  one  moves  downstream.  Typical  such  plateaus  are  seen 
in  cross-sectional  views  of  Figure  20  shown  in  Figure  - 1 • 

The  calculations  outlined  above  demonstrate  that,  although  quantitatively 
different  for  the  axisymmetric  hypersonic  wake,  an  essential  feature  of  the 
wake-at-incidence  problem  could  indeed  be  obtained  in  our  wind-tunnel.  As 
described  in  the  following  paragraphs,  the  shear-layer  feature  was  clearly 
evident  from  the  data  collected  ao  far. 


A. 3.2  CHOICE  OF  INCIDENCE  SETTINGS 


To  study  the  effect  of  incidence,  more  than  one  incidence  angle  should  be 
chosen,  if  possible;  a  judicious  choice  of  angles  is  needed  in  order  to 
illuminate  the  effect  without  undue  labor.  It  is  known  on  theoretical  grounds 
that  beyond  the  detachment  angle  (3A  degrees)  the  inviscid  wake  freezes  to 
a  configuration  not  very  sensitive  to  angle  of  incidence.  lests  verifie  ^ 
this  effect  and  i  the  detached  case  it  was  decided  to  work  only  with  a  -  U 
(plate  normal  to  the  stream)  so  that  the  blunt-body  turbulent  wake  could  at 

least  be  probed. 


In  the  pre-detachment  range  of  a  wake  "pitot  maps",  suchas  shown  for 
a  *  0°  on  Figure  30  of  Reference  1,  were  taken  at  a  =  10  ,  15  ,  20  and  JU 
Such  maps  are  easy  to  construct  and  are  invaluable  in  giving  an  overall 
picture  of  the  flow  field.  Since  the  shear-layer-like  behavior  was  present 
in  each  case,  the  angles  d=  10°  and  20°  were  chosen  for  further  study. 
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FIGURE  18  PHOTOGRAPH  OF  WIND-TUNNEL  TEST-SECTION  SHOWING 
TWO-DIMENSIONAL  MODEL  SUPPORT 
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FIGURE  19  SCHEMATIC  OF  TWO-DIMENSIONAL  WAKE  NOMENCLATURE 


In  designing  the  experiment,  the  farthest  downstream  distance  where  measure¬ 
ments  can  be  made  is  important  to  know.  In  the  wind-tunnel  this  distance  is 
shortened  when  the  model  is  pitched  because  of  the  earlier  reflection  of  the 
strengthened  shocks  from  the  tunnel  wall.  In  the  present  set-up  the  farthest 
downstream  distance  fora  -  0° ,  10°,  20°  and  90°  is  7.8  inches,  7  inches, 

6.2  inches  and  <4.7  inches,  respectively.  Note  that  for  a  ■  0°,  a  distance 
of  7.8  inches  implies  a  distance  of  about  1,800  virtual  body  diameters. 

A.4  DIAGNOSTIC  INSTRUMENTATION 

ihe  diagnostic  tools  used  for  this  experiment  were  identical  to  those  employed 
for  the  WED  (see  Reference  1).  The  0.00005  inch  diameter  Pt  10  percent  Rh 
hot-wire  used  for  this  work  managed  to  stay  Intact  through  the  rean-flow  and 
turbulence  measurements  and  through  most  of  the  spectral  measurements.  This 
remarkable  wire  survived  45  calendar  days  of  use,  during  which  time  it  was 
operated  for  about  70  hours  in  supersonic  flow.  Through  frequent  oven-  and 
flow-calibrations  small  changes  in  its  characteristics  could  be  monitored  and 
accounted  for. 

4 . 5  PROCEDURE 

This  experiment  was  designated  WEF  (Wake  Experiment  F)  and  the  a  *  10°,  20°, 
and  90°  sub-experiments  were  designated  WEF-10,  WFF-20,  and  WEF-9Q,  respec¬ 
tively.  As  done  with  the  WED  experiment,  the  distance  between  the  model 
mid-chord  and  the  end  (farthest  downstream)  point  was  marked  by  X-stations 
spaced  by  0,300  inch  increments.  Continuous  analog  traverses  were  taken  at 
each  X  station  with  the  pitot  and  static  probes  and  the  hot-wire  anemometer, 
the  latter  operated  either  in  the  mean-flow  or  the  turbulent  mode.  From 
these  data  the  mean  and  turbulent  properties  can  be  extracted. 

4.6  RESU1.TS  OF  MEASUREMENTS 

As  already  mentioned,  the  actual  viscous  wake  measurements  were  preceded  by 
tracing  the  inviscid  wake  in  order  to  observe  the  overall  flow  field  and  to 
examine  the  existence  of  the  "shear-layer"  effect.  Figure  22  compares  a 
typical  pitot  traverse  with  the  pitot  pressure  expected  from  the  character¬ 
istics  computation  shown  on  Figures  20  and  21  fora  =  20°.  It  is  obvious 
from  the  measurement  that  two  shocks  appear  on  each  side  of  the  inviscid 
wake,  instead  of  one  as  given  by  the  characteristics  method.  On  the  wind¬ 
ward  side,  and  in  addition  to  the  bow  shock  at  the  leading  edge,  a  second 
shock  is  thought  to  originate  in  the  recompression  region  of  the  wake  "neck", 
as  illustrated  tentatively  on  Figure  19.  In  the  leeward  side  a  shock  appears, 
of  course,  at  the  model  leading  edge  in  addition  to  the  expansion  fan.* 


*Photography  of  the  two-dimensional  model  itself  is  impossible  because  of 
obstructions.  The  ident  i  f  icat ion  of  the  additional  two  shocks  is  left  to 
reasonable  conjecture. 
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Nevertheless  the  agreement  shown  on  Figure  21  between  the  calculated  and  ^ 

measured  pitot  pressures  is  adequate;  both  demonstrate,  in  fact,  the  plateaus 
on  either  side  of  the  turbulent  wake. 


A  second  preliminary  measurement  made  before  detailed  studies  were  attempted 
concerned  the  existence  of  turbulence.  Qualitative  traverses  with  the  hot¬ 
wire  showed  that  angle  of  attack  did  not  inhibit  the  appearance  of  turbulence 
in  the  inner  wake.  Proceeding  to  the  latter,  Figure  21  shows  pitot-pressure 
traces  at  typical  X  stations  for  all  three  experiments,  i.e.,  WEF-10,  WEF-20, 
and  WEF-90.  A  trace  for  WED  (corresponding  to  WEF-0)  is  also  shown.  The 
object  of  this  figure  is  to  demonstrate  the  wake  symmetry  which  understandably 
obtains  for  the  symmetric  body  configurations,  namely  for  WED  and  WEF-90. 

Note,  also,  the  expected  change  in  "f ree-stream1  pitot  level  tor  the  WFF-10 
and  WEF-20  traces,  illustrating  the  "shear-layer"  effect. 


Typical  static  pressure  traverses  are  shown  on  figure  24.  Although  these 
traces  appear  reasonably  smooth,  it  should  be  recalled  that  the  static-probe 
resolution  is  quite  poor,  since  the  diameter  of  the  probe  was  never  smaller 
than  five  times  or  so  the  wake  thickness.  On  the  other  hand  there  is  strong 
belief  that  the  mixing  in  the  far  wake  occurs  at  constant  pressure,  or  else 
the  whole  flow  system  should  shift  progressively  in  the  directly  of  decreasing 
pressure.  In  fact,  this  probe  gave  seemingly  good  resolution  in  regions  of 
sensible  pressure  gradients,  such  as  at  or  near  shock  waves. 


Hot-wire  measurements  consist  oi  the  usual  four  sets:  one  set  onsists  ■  i 
five  curves  at  each  X  station,  necessary  to  complete  the  inputs  for  the 
mean  flow  measurements.  A  second  set,  with  15  curves  (each  at  different 
current)  traces  the  variation  of  mean  voltage  across  the  wake;  a  companion 
set  of  15  curves  plots  the  mean-square  of  the  wire  a.c.  'turbulent)  component. 
Here  one  should  note  that  the  turbulent  fluctuations  are  highest  on  the 
side  of  the  wake  where  the  velocity  gradients  are  the  highest.  Spectral 
measurements  have  not  been  taken  as  yet. 
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FIGURE  24  TYPICAL  STATIC-PROBE  TRACES  OF  THE  TWO-DIMENSIONAL  WAKE  FOR  A  20-DEGREE  INCIDENCE  ANGLE 
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SECTION  V 

SUPERSONIC  WAKE  WITH  HE/T  TRANSFER 


5 . 1  INTRODUCTION 

There  are  several  distinct  reasons  why  the  pressure  of  heat  transfer  is 
important  to  understanding  the  fluid  mechanics  oi  the  high-spec  wa  c 
In  the  first  place, it  is  well  known  from  stability  arguments  (Reference  14) 
that  the  location  of  wake  transition  to  turbulence  is  sensitive  to  heat 
exchanged  between  the  body  and  the  flow.  Secondly,  the  extent  to  which 
heat  transfer  changes  the  self-preservation  behavior  is  not  lG^  transfer- 

Third  the  fluctuations  in  the  wake  are  much  more  complex  with  heat  transfe  , 
for  example,  the  sample  relationship  between  velocity  and  density  fluctua¬ 
tions  valid  for  the  isoenergetic  wake  does  not  apply. 

In  the  present  series  of  experiments,  it  is  planned  to  study  the  mean 
intermittent  and  turbulent  flow  component  behind  an  axi-symmetnc  coo  ed 
model  and  a  two-dimensional  heated  model.  The  latter  experiment  has  been 
initiated  and  most  of  the  test  have  actually  been  completed.  It  will  be 
described  in  this  section. 

5.2  MODEL  AND  FACILITIES 

This  experiment  is  performed  in  the  same  wind-tunnel,  at  the  ‘same  condi tion, . 
with  the  same  model  and  concurrently  as  the  anglc-of-at each  experiment 
discussed  in  Section  IV.  Provisions  to  the  model  support  additional 
described  in  the  latter  section  as  neressory  for  the  angle-of-attach 
measurement,  were  incorporated  specifically  for  the  heat- transfer  me  su 

ment . 

The  steel  ribbon  is  hentod  electrically  by  a  Harrison  Laboratories 
Model  No.  81AA  regulated  dc  power  supply.  Load  regulation  is  c  to  c  o 
0.05  percent, and  power- line  ripple,  already  low,  is  node  o  cc  vc  y 
imperceptible  because  of  high  tncrmal  inertia  of  the  boated  ribbon.  The 
ribbon  current  con  be  adjusted  from  0  to  25  amperes  dc  and  is  continuously 
displayed  on  a  Spoedomax-W  strip-chart  recorder  running  at  a  speed  of 
1  inch  of  paper  per  second.  Aided  by  the  chart  resolution  and  a  fine 
variable  resistor  added  externally  to  the  circuit  the  operator  can  thus 
readjust  and  control  the  ribbon  current  to  about  0.015  percent  of  the 
desired  value,  an  error  lying  much  below  the  current  change  needed  for 
perceptible  changes  in  the  heated  wake  characteristics.  The  power  settled 
upon  for  this  measurement  will  be  given  farther  below. 

Expansion  of  the  model  by  thermal  relaxation  necessitated  a  scheme  by 
which  the  ribbon  torsion  could  be  adjusted  automatically.  This  was  done 
by  spring- loading  one  end  of  the  ribbon  so  that  the  model  could  expand 
when  heated  and  retract  to  a  stretched  position  when  the  electrical 
power  was  turned  off. 
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5.3  ADJUSTMENT  OF  OPERATING  CONDITIONS 

To  bring  about  easily  observable  heat-transfer  effects,  it  was  desired  to 
dissipate  the  highest  possible  amount  of  power  in  the  ribbon.  Several 
othei  conditions  should  be  satisfied  simultaneously;  the  ribbon  should  not 
weaken  structurallv  or  deform,  the  turbulence  should  not  be  inhibited 
the  poiniwhere  the  wake  became  wholly  laminar,  and  the  ribbon  temperaure 
non-uniformity  because  of  heat  conduction  into  the  supports  Should  be  kept 
minimal .  There  was  furthermore  some  fear  that  the  increased  wa  e  emp 
tures  could  restrict  the  operating  temperature  of  the  hot-wire  or  even 
destroy  it  when  both  the  model  and  the  hot-wire  operating  temperatures 
were  high.  As  it  developed,  all  these  criteria  were  satisfied  once  the 

ribbon  maximum  power  was  set. 

Tii .» c  was  accomplished  by  gradually  increasing  the  ribbon  current  until  Lhc 
mo  I  g  orC™  ly  in  the  supersonic  How.  The  current  was  then  reduced 
"til  the  glow  disappeared.  At  this  point  the  ribbon  “treat  was  set  at 
19.90  amperes  giving  a  voltage  drop  across  i to  4; 344  vo  s  „ 

rntnl  nower  dissipation  ol  W  -  oo.o  wants.  XL  "  .  .  \ 

Ih"  pLwer  to  the' "Kinetic"  power  of  the  flow 

Utilizing  the  viscous  drag  coeftictent  of  the  ribbon  from  the  • 
incidence  WED  experiment  (Section  II)  a  ratio  ot  W/P  -  1.  i»  obtainc 

After  the  necessary  tests  for  finding  the  true  zero-twist  and  zero-incidence 
angles  (sec  Reference  1),  the  hot-wire  was  employed  qualitatively  t 
studying  the  transition  to  turbulence  with  and  the  transverse  (span-w  se) 
uniformity  of  the  heating  process.  Traasicioa  was  detected  by  the  usual 
method  of  noting  the  poinl  of  maximum  turbulence  intensity  along  the 
center  Une  ^ferencc'  15).  Figure  25  shows  the  remarkable  dl  ferences  in 
transition  distance  observed  with  heating*  this  distance,  located  abo 
100  inches  virtual  diameters  (momentum  thicknesses)  behind  the  adiabatic 
body  moved  to  about  800  diameters  when  heating  was  applied.  This 

ith,havallalle£«ngeCiiA)the  axial'diwcti^o™"^*^  dlagnoslics  could  be 

i  .  r  hnlf  the  ranee  available  for  wake  studies  with  the 

For  this  reason  probings  began  again,  as  in  wuu,  3 

body  rather  than  at  the  transition  distance. 

Visual  observation  of  the  ribbon  Lollowing  its  overheating  to  the  glowing 
no in showed V that  its  rear  half  (from  about  mid-chord  to  the  trailing 
edge)  had  been  discolored  in  the  way  typical  of  metals  alter  exposure 
higl  temperatures.  That  the  ribbon  temperature  increases  -long  it ^  ^ord 
toward  the  trailing  edge  is  fully  expected  because  of  he  more  efficient 

^Henceforth"  "heating"  refers  to  the  86.8  watts  settled  upon  as 
noted  previously. 
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TRANSITION  MOVED 
BY  HEATING 


FIGURE  25  DELAY  OF  TRANSITION  TO  TURBULENCE  IN  THE  WAKE  BROUGHT  ABOUT  BY  HEAT  TRANSFER 


heat  transfer  for  the  thinner  boundary  layer.  The  discolored  region 
extended  uniformly  along  the  span,  but  it  was  obseived  to  stop  abruptly 
about  one  ribbon  chord  before  it  reached  the  tunnel  sidewalls.  The  ends 
of  the  ribbon  were  therefore  cooler  since,  of  course,  heat  conduction  to  the 
solid  supports  increased  at  those  points. 

The  observations  described  above  are  crucial  to  the  span-wise  uniformity  of 
heat  transfer  from  the  ribbon  to  its  wake,  that  is,  tne  question  of 
true  two-dimensionality  of  the  experiment.  This  question  was  further 
examined  by  taking  hot-wire  profile  of  the  local  total  wake  temperature  at 
fixed  distances  behind  the  model  and  different  spanwise  stations.  Within 
an  inch  of  the  spanwise  direction  around  the  center-span  position 
(covering  about  35  percent  of  the  total  span) ,  the  wake  axis  temperature 
varied  by  about  10  percent;  at  least  part  of  this  variation  is  attributable 
to  slight  flow  non-uniformities.  The  two-dimensionality  was  thus  judged 
to  be  adequate. 

5.4  PROCEDURES  AND  INSTRUMENTATION 

The  instrumentation  used  for  this  work  was  identical  to  that  used  for  the 
WED  and  WEF  experiments.  Although  the  ribbon  itself  was  heated  highly, 
turbulent  diffusion  reduced  the  local  stagnation  temperatures  to  the  point 
where  no  question  of  damage  to  the  probes  by  high  temperatures  ever  arose. 

In  fact  the  0.00  05  inch  diameter  wire  was  operated  in  the  same  full  range 
of  heating  currents  (maximum  of  about  7.5  mi 1 liamperes)  without  ever 
failing  by  causes  attributable  to  high  wake  temperatures. 

The  present  experiment,  done  with  the  heated  two-dimensional  body  at  zero 
incidence,  is  designated  WEG  (Wake  Experiment  G)  .  The  measurements  were 
performed  at  each  of  23  axial  positions  along  the  wake.  (XSTATIONS  1 
through  23)  beginning  at  0.658  inch  from  the  ribbon  mid-chord  point.  At 
each  station  one  pitot-pressure  transverse,  one  static-probe  transverse, 
and  five  hot-wire  traverses  were  performed  to  collect  data  sufficient  to 
describe  fully  the  mean  flow  field.  Each  wire  traverse  was  of  course 
performed  at  different  heating  currents.  For  the  turbulence,  15  traces  of 
mean  wire  voltage  and  rms  voltage  were  collected  at  each  station.  For  the 
spectra,  ten  distinct  radial  locations  were  chosen  at  which  spectra  were 
recorded  at  each  of  six  different  heating  currents.  Since  the  acquisition 
and  reduction  of  the  spectral  data  is  extremely  time-consuming,  data  were 
taken  only  at  alternate  XSTATIONS  (the  odd-numbered  stations  1,  3,  5,... 23). 
Thus  a  total  of  10  x  6  x  12  =  720  spectra  were  recorded. 

It  is  perhaps  advisable  to  reiterate  here  the  data  reduction  process 
through  which  the  WEG  data  are  now  progressing.  The  mean- flow  data  are 
digitized  and  put  into  the  WEB-V  and  WEB-VIII  computer  programs  which 
have  for  some  time  been  prepared  for  the  Philco-2000  digital  computer. 
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The  spectral  data  are  first  used  into  the  WEB-II  program  from  which  one 
extracts  the  so-called  error  ratios,  that  is,  a  correction  factor  necessary 
to  use  in  correcting  the  system  frequency  response.  The  output  of  this 
program  and  the  turbulence  data  are  then  used  into  the  WEB-III  and  WEB-IV 
programs  to  deduce  the  turbulence  properties  and  further  into  the  WEB-VII 
program  from  which  turbulence  properties  in  scaled  form  are  obtained. 

5.5  RESULTS  TO  DATE 

The  raw  data  obtained  to  date  give  the  following  indications  of  wake 
behavior : 

(1)  Profiles  are  symmetric  as  expected.  This  applies 
equally  to  the  pilot,  static  and  hot-wire  data. 

(2)  As  expected  when  heat  transfer  is  at  play,  the 
flow  profiles  at  the  far  wake  positions  stand 
out  much  more  clearly  than  in  the  adiabatic 
experiment.  The  local  total  temperature,  especi¬ 
ally,  barely  distinguishable  from  the  free-stream 
value  in  WED,  has  an  outstanding  intensity. 

(3)  The  turbulence  intensity  increases  at  anv  given 
station  following  transition. 

(4)  The  spectra  again  shew  the  energy  prominence 
they  also  displayed  in  WED;  the  peak,  however,  is 
displaced  to  lower  frequencies,  perhaps  reflecting 
an  overall  thickening  of  the  wake. 
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APPENDIX  I 

DATA- HANDLING  PROCEDURES 


1.1  INTRODUCTION 

In  the  work  done  on  the  turbulent  wake  structure  by  this  company  in  the 
recent  past,  it  lias  been  demonstrated  that  automated  procedures  in  die  wind- 
tunnel  have  enabled  a  very  great  volume  ol  relevant  data  to  be  acquired. 

It  is  necessary  to  accelerate  the  reduction  rate  of  these  data  between  the 
point  of  transducer  output  and  the  point  of  input  into  the  computer. 
Furthermore,  it  is  clear  that  the  present  hot-wire  anemometer  system  ought 
to  be  improved  to  handle  higher  frequency  ranges.  Finally,  it  is  necessary 
to  re-direct  the  experimental  equipment  for  near-future  requirements  in¬ 
volving  novel  measurements  ol  turbulence.  Considerable  effort  lias  been  put 
into  the  planning  of  new  data-produc ing  and  handling  techniques. 

1.2  DATA  ACQUISITION:  HOT-WIRE  ANEMOMETRY 

As  Paragraph  2.4  explains,  the  present  hot-wire  anemometer  system  is  capable 
of  yielding  data  to  frequencies  of  order  1MHz.  A  new  system  is  being  planned 
to  extend  this  range  to  several  MHz.  This  system  is  based  on  the  use  of  an 
extremely  fine  hot-wire,  0.00001  inches  in  diameter,  and  about  1  .''0.  inches 
long.  The  time-constant  of  this  wire,  unaided  by  electronic  compensation,  will 
be  between  10  and  20  microseconds.  The  wire  signal  will  be  directed  into 
a  new  constant -cur rent  hot-wire  amplifier  which  has  now  been  designed.  This 
amplifier  will  utilize  modest  amplification  ratios  and  a  f loor-Lo-ce i 1 ing 
ratio  not  much  in  excess  of  500.  Solid-state  amplifiers  in  several  stages 
will  produce  a  flat  s tra ight-ampl i f ier  response  to  beyond  1  MHz  and  a  com¬ 
pensation  capability  in  excess  of  that  frequency.  Self-contained  mercury-cell 
power  supplies  will  heat  the  wire.  Two  such  amplifiers  are  being  planned. 

The  next  element  in  the  line  will  be  a  high- frequency-response  (2  MHz  or 
better)  magnetic  tape  recorder  with  a  minimum  of  seven  channels,  any  number 
of  which  will  be  preselected  to  accept  dc  signals,  with  the  remainder  used 
for  the  high-frequency  ac  signals.  The  recorder  will  be  used  for  a  variety 
of  important  functions.  In  a  typical  wake  experiment  two  ac  channels  will 
carry  the  direct  hot-wire  output  (the  second  channel  being  active,  for  ex¬ 
ample,  in  correlation  measurements  with  two  wires),  two  dc  channels  will  be 
carrying  voltages  proportional  to  two  position  coordinates,  and  the  re¬ 
maining  three  dc  channels  the  wire  rms  output,  wire  mean  voltage  and  wire 
current.  By  recording  at  closely-spaced  points  and  different  wire  currents, 
it  is  hoped  to  obtain  a  complete  record  of  an  entire  experiment  stored  in 
two  15  inch  reels  of  magnetic  tape. 

The  signals  most  time  consuming  operation  in  the  current  work,  both  as 
regards  the  data  acquisition  and  its  digitizati  n,  is  the  measurement  of 
spectra.  Currently,  several  hundred  spectra  are  obtained  for  each  experiment. 
Each  spectrum  is  obtained  in  analog,  bivariate  form  and  then  digitized  by 
hand  for  key-punching  into  computer  card  format.  It  is  planned  to  bypass 
this  operation  by  taking  the  ac  signal  from  the  magnetic  tape,  digitizing 
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it  electromechanically  and  computing  its  spectrum  by  using  a  computer  pro¬ 
gram  (the  "Fast  Fourier  Transform"  program). 

The  remainder  of  the  operations  needed  tor  mean  and  turbulence  measurements 
involve  only  dc  measurements  e.g.,  a  dc  transducer  output  voltage  and  a 
dc  position  voltage  which  can  be  obtained  by  either  of  three  methods: 

(1)  from  the  points  of  the  magnet ic -tape  file  of  the  complete  experiment 

(2)  by  a  second  file  easily  made  by  repeating  the  experiment  by  continuous 
traverses,  or  (3)  by  a  method  identical  to  (2)  but  where  the  two  voltages 
are  digitized  and  punched  (on  paper  tape  or  cards)  in  real  time  without 
recourse  to  magnetic  tape.  In  either  case  large  savings  in  time  ar 
possible . 
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14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report  Key  words  must  be 
selected  so  that  no  security  classification  is  required  Identi¬ 
fiers.  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional 
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